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We consider a superconducting quantum point contact in a circuit quantum electrodynamics setup. We study 
three different configurations, attainable with current technology, where a quantum point contact is coupled 
galvanically to a coplanar waveguide resonator. Furthermore, we demonstrate that the strong and ultrastrong 
coupling regimes can be achieved with realistic parameters, allowing the coherent exchange between a super- 
conducting quantum point contact and a quantized intracavity field. 



Introduction. — A quantum point contact (QPC) is a super- 
conducting constriction, where electronic transport is sup- 
ported by small number of conducting modes 1 9 , and the size 
of the constriction is small compared to the inelastic mean 
free path LLM<il Above the superconducting critical tempera- 
ture, the normal state resistance of the contact is determined 
by the transmission of those channels. In the superconduct- 
ing state, the dissipationless current is carried by the Andreev 
bound states (ABS) localized at the junction region over the 
superconducting coherence length^. There is a pair of ABS 
for each conducting channel and their energies, ±£U depends 
on the phase difference cf> across the junction and the trans- 
mission coefficient r of the conducting channel Ea(<I>, t) = 



Ay 1 — r sin 2 (0/2), where A is the order parameter of the 
bulk superconducting electrodes. The energies of ABS are 
defined with respect to the Fermi level, and the current car- 
ried by each state is given by I± = ±(2e/ti)dEA/dcf). The 
total dc Josephson current through a QPC, consisting of a sin- 
gle conducting channel, is given by the contribution of each 
ABS weighted by the state occupation, / = I-ti(—Ea) + 
I+ii(Ea), where n(E) is the Fermi distribution function. Re- 
cent experiments on QPC have shown the correctness of de- 
scribing the Josephson effect in terms of the ABS™. In ad- 
dition, a direct spectroscopy of the ABS was done in a carbon 
nanotube attached to superconducting electrodes by means of 
tunneling spectroscopy^. 

On the other hand, several theoretical works have studied 
the ABS and their dynamics^" 9 -! In particular, it was recently 
proposed to use a QPC as a qubit 4 , introducing an effective 
Hamiltonian describing the two-level system of ABS. It has 
also been shown that in the zero-temperature limit the dynam- 
ics of the QPC in a microwave field is well described by a two- 
level effective Hamiltonian 9 . In this sense, it is predicted that 
the transition between the levels induced by the microwave 
field leads to a drastic suppression of the dc Josephson cur- 
rent. However, though several attempts to detect such transi- 
tions have been carried out, these have not been directly ob- 
served yet. In this context, it would be of great interest to 
develop a theory involving a QPC coupled to a superconduct- 
ing resonator in the context of state-of-the-art circuit quantum 
electrodynamics (QED) 18 '' 21 [ This may serve as an alternative 



physical scenario for probing the ABS dynamics, as well as a 
testbed for fundamental aspects in light-matter interaction. 

In this work, we propose a circuit QED setup using a super- 
conducting QPC coupled galvanically 22 ^ to a cop lanar waveg- 
uide resonator (CWR). Following previous works^GS we de- 
velop a theoretical framework to study three experimental 
configurations shown in Fig.[T| and we present numerical sim- 
ulations showing distinct cavity QED features as compared to 
the semiclassical case. Furthermore, we demonstrate that, for 
the configuration in Fig. lc, one can reach the ultrastrong 
coupling (USC) regime of circuit QED^HSZl 

The Model. — Our generic setup consists of a superconduct- 
ing loop containi ng a QP C coupled galvanically to a inhomo- 
geneous resonato r 22 ! 26 ! 27 ! (see Fig. fla). We neglect the geo- 
metrical inductance of the loop, which we consider 15 of the 
order of 20 pH, with respect to the inductance of the QPC 
(Lqpc — 10 nH^. The system Hamiltonian reads 



H = H, 



QPC I 



H, 



CWR- 



(1) 



We assume cryogenic temperatures and external frequencies 
smaller than A — Ea- In this case, one can use the effective 
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FIG. 1. (Color online) (a) Superconducting loop containing a 
QPC coupled galvanically to a CWR (see inset), (b) Atomic 
SQUID (aSQUID) with a QPC in series with a Josephson junction, 
(c) aSQUID with the QPC in parallel with a Josephson junction. 
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Hamiltonian derived by Zazunov et al^to describe the QPC, 

^QPc(0a) = Ae~^^ /2 [cos(0 Q /2)a z + Vfsin(0 a /2)a,] 

(2) 

where r = 1 — r stands for the reflection coefficient and the 
Pauli matrices a x y z are written in the ballistic basis, defined 
by the eigenstates of the current for a perfectly transmitted 
conducting channel. 

The CWR is described by a sum of harmonic oscillators 
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such that the total phase distribution on the CWR reads 
4>{x,t) = ~Yl, n u n {x)(j) n {t), where u n (x) is the spatial eigen- 
mode. The phase variation along the segment of the res- 
onator Ax, shared with the QPC loop, see Fig. [T^, reads 
4> r = J2 n 3n<t>n, where S n stands for the variation of the 
nth spatial eigenmode, u n {x), along Ax. The phase (f> n 
is related to the annihilation and creation operators <j> n — 
(2e/K)wh/ {2ui n C r ){a\ l + a n ), where uj n is the correspond- 
ing eigenfrequency, and C r is the capacitance of the CWR; the 
conjugate momentum reads 9 n = i(h/2e)y/h(jj n C r /2(a'l l — 
a n ). Typical resonator parameters are C r ~ 850 fF, the 
impedance Z = y/L r /C r ~ 50 fl, where L r is the CWR 
inductance, and frequency u> r /2ii ~ 1 — 10 GHz for the fun- 
damental mode^. 

In all cases displayed in Fig. [T^, the description of the 
CWR in terms of eigenmodes is still valid if the inductance 
L x of the resonator segment Ax is small as compared to 
the inductance of the QPC. Otherwise, the mode structure 
will depend strongly on the nonlinearities coming from the 
QPC current. Under this approximation most of the current 
will flow through the resonator, and the QPC will act as a 
small perturbation. For a CWR made of aluminum with width 
S = 50 nm (at the constriction), thickness of the central elec- 
trode t = 50 nm, length of the constriction Ax = 5 urn and 
distance between the ground plane and the edge of the central 
line W — 4.95 /im, we estimate L x ~ 8 pH <C £qpc> thus 
justifying the approximation. In order to simplify our anal- 
ysis, we consider the simplest case where the QPC interacts 
with the lowest frequency eigenmode, n — 0, supported by 
the resonator. 

The CWR and QPC are mutually coupled via the phase 
difference along the segment of the resonator Ax. Due to 
the flux quantization in superconducting loops, one can re- 
late the external flux threading the loop, <£>, with the super- 
conducting phase differences, fn — 27r/ + 2ttN, where 
we have defined the frustration parameter / = being 
$o = h/2e the flux quantum. In particular, for the situa- 
tion of Fig. [T^, the phase across the QPC (<p a ) and the phase 
drop along the segment Ax are related by 4> a = <fr r + 2vr/. 
The superconducting phase difference <f> r is given by <f> r — 
S (2e/H)(H/2ui C r ) 1 ^{a + at ). We will also assume that 
\4> r \ <C 1, a condition satisfied with current available val- 
ueJSl. In this case, we can expand the Hamiltonian §2§ up 
to first order on <j> r . In addition, we rotate Hqpc m Eq. (|2| 
to the instantaneous Andreev states basis? by means of the 
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FIG. 2. (Color online) (a) Contour plot showing population inversion 
of Andreev levels for initial state j^o) = |— }jl)r- The inset shows 
the mean value of <j z as a function of normalized time, t(l/go), at 
the resonance condition fo = 0.4712. 



transformation H, 



(A) ( 
QPCV 



„) = UH^cttaWK where U = 



QPC 



e i<Pa<? y e i 4 Under this unitary transformation the whole 
Hamiltonian reads 

H = hu> Q a) a+E A (f,T)a z -hg(a+a' !< )(A z a z -A x a x ), (4) 



where g 



1 Al 

2 h 



im(TTf)5 (2e/K)(h/2u} C r ) 1 / 2 is the cou- 



pling strength, while A z and A x are the longitudinal and 
transversal coefficients, respectively, 



Az - E A (f, t)' K 



y/fsin(7r/) 
E A (f,r) ' 



(5) 



Notice that for a frustration / ~ 0.5, a large (small) contri- 
bution of the transversal (longitudinal) coupling is obtained. 
This is important since in a small vicinity of / = 0.5 the QPC 
and the resonator can exchange excitations, so that ABS spec- 
troscopy is possible. On the other hand, away from this work- 
ing point the QPC -resonator interaction is dominated by the 
longitudinal coupling (A z ), where the field of the resonator 
will be displaced depending on the state of the QPC. These 
features show the versatility of our setup in order to engineer 
the QPC-resonator coupling, and will be the base for studying 
the underlying physics in three different setups. 

QPC coupled to a coplanar waveguide resonator. — We first 
analyze the configuration of Fig. [T^, i.e. a superconduct- 
ing loop containing a QPC galvanically coupled to the CWR. 
The dynamics of this setup is described by the Hamiltonian 
in Eq. Q. This resembles the coupling of a spin to a mag- 
netic field with axial and longitudinal components and has the 
same structure like th e Ham iltonian describing Cooper pair 
boxes or flux qubits^- 24 * 26 l 27 l It is worth mentioning that for 
parameters satisfying the conditions: ojq + 2Ea(J,t)/H 3> 
{gA x , \u>q — 2Ea(J ,r)/h\} and wo 3> gA z , the Hamilto- 
nian (|4]i can be approximated by the the Jaynes-Cummings 
mode™ for which the rotating wave approximation holds. 
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We have studied numerically the dynamics of the Hamil- 
tonian in Eq. Q for an initial state having a single excita- 
tion in the resonator and the QPC in the ground state, i.e. 
|^o) = | — ) |l) r - Figure [2] shows a contour plot of (a z ) as 
a function of a normalized time, and the external flux ap- 
plied to the loop close to the resonance condition, ojq = 
2(A/ft)(l - rsin 2 (7r/)) 1 / 2 , that leads to the working point 



fo = — arcsin 

7T 



V 2A / 



1/2- 



(6) 



This figure reveals the population inversion of Andreev lev- 
els near the working point fo, as expected from the Jaynes- 
Cummings dynamics. This feature is more clearly shown 
in the inset of Fig. |2j where we show the time evolution 
of (<j z ) at the resonance fo = 0.4712. In this simula- 
tion the CWR is described by a single mode with frequency 
cj /27r = 10.52 GHz, and a phase drop \<f> r \ ~ 0.0013 
(see Ref. l23l . For the QPC made of aluminum, we con- 
sider only one conducting channel with high transmission, 
t = 0.994, and A/ft ~ 44.256 GHz (corresponding to a 
gap A ~ 0.183 meV for Alj^. These values lead to a ratio 
Sto/wo ~ 0.002, where go = ^TS (2e/h)(h/2Lo C r ) 1 / 2 . 
In principle, for these particular values the rotating-wave ap- 
proximation can be safely used. 

It is worth mentioning that the vacuum Rabi oscillations 
shown in Fig. [2] involves a time domain measurement of the 
Andreev levels. This can be done by, firstly, applying an ex- 
ternal field via an on-chip flux line on the loop containing 
the QPC to far detune it from the resonator. Then, an am- 
plitude shaped flux pulse /d c is applied to tune the qubit into 
resonance with the cavity field for a variable time r. Then, 
probe the cavity frequency in a dispersive regime where the 
shift of the resonator frequency is proportional to the Andreev 
level population imbalance 29 . In addition, the spectrum of 
the QPC-resonator system is shown in Fig. [3] that should be 
obtained via a cavity transmission measurement at frequency 
The inset shows the avoided crossing between the sin- 
gle photon resonator state and the excited Andreev level as the 
QPC is brought into resonance at fo = 0.4712. 

Atomic SQUID coupled to a coplanar waveguide 
resonator. — For a proper characterization of the QPC 
and to control its phase, experiments have been carried out on 
an asymmetric SQUID loop consist ing of a Josephson tunnel 
junction (J J) in series with the QPCPEEH. Here, we consider 
such an atomic SQUID (aSQUID) galvanically coupled to 
the resonator as shown in Fig. [IJ>. In this configuration, the 
relation among the phases associated to the QPC (4> a ), the 
Josephson junction ((j> p ), and </> r reads <j ) a = 4>r + 4>p + 2vr/. 
The experiments considering this setup 12 13 15 , implement a 
large Josephson junction such that the Josephson energy is 
larger than the charging energy, Ej Ec- In this case, 
the junction behaves as a harmonic oscillator that will be 
described by annihilation and creation operators b,b\ and 
plasma frequency lo p = (2elc /ftCj) 1 / 2 with Cj and Iq 
the Josephson junction capacitance and critical current, 
respectively. Within this description the Hamiltonian for the 
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FIG. 3. (Color online) Energy spectrum for the configuration showed 
in Fig. [T^. The color codes stand for the energy differences — 
computed numerically — blue-line stands for (E\ — Eq)/H, green- 
line for (E2 — Eo)/k, red-line for (E3 ~ Eo)/h, and light-blue-line 
for (E4 — Eo)/h. The inset shows the characteristic anticrossings 
representative of the Jaynes-Cummings model. 



setup of Fig.[TJ) reads 

H = hujQ(ra + huj p b^b + Ea{J, t)ct z 
A 2 

- t— sm(irf)(A z <r z - A x (j x )(4> r + <f> p ) 



A 2 



2y/r sin(7r/) cos(7r/)cr x 



(cos 2 (tt/) -rsin 2 (7r/))cr z 



(7) 



where we have defined <j) p = {2e/h)(h/2C J uj p ) 1 / 2 {b + 

In obtaining the Hamiltonian Q, we considered realistic 
values of <fi r and <f> p up to second order. For instance, by 
taking a Josephson junction with critical current Ic ~ 1/iA, 
capacitance Cj ~ 100 fF, and plasma frequency u> p /2tt ~ 
27 GHz 30 , one obtains \<f> p \ ~ 0.16, while for the the phase 
drop at the constriction we take \cf) r \ ~ 0.0013. Notice that 
the dynamics results into a Stark shift of the Andreev lev- 
els depending on the number of excitations in the JJ and the 
CWR. Figures and show the spectrum for the setup of 
Fig.fTJ), calculated via a numerical diagonalization of Hamil- 
tonian ([7]). It shows that the presence of the JJ induces an 
ac-Stark shift of the Andreev levels as depicted in Fig. |4j). We 
have estimated numerically a renormalized resonance condi- 
tion, fo ~ 0.47059 where the QPC exchanges excitations with 
the resonator, as observed in the displayed anticrossings. 

Ultrastrong coupling regime. — Let us consider the setup 
shown in Fig. [T]:, where the Josephson junction is now part 
of the resonator line. This setup is characterized with two 
main effects: (i) the renormalization of resonator eige nfre- 
quencies, which are now given by the dispersion 

relatiotPSl 

k n = (2Lo/Lj)(l — cj 2 /w 2 ) cot(k n l), where Lj is the junc- 
tion inductance, Lq the inductance per unit length of the res- 
onator, and I is half of cavity length, (ii) An increase of the 
coupling strength g due to the local modification of the indue- 
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FIG. 4. (Color online) Energy spectrum for two cases of Fig.^ (a) 
and (b) correspond to Fig. lb, while (c) and (d) to Fig. lc. The 
color codes represent the computed energy differences. Blue-line 
stands for (Ex — Eo)/H, green-line for (E2 — Eo)/h, red-line for 
(E 3 - E )/h, and light-blue-line for (E 4 - E )/H. 



tance of the single-mode resonator. For instance, if one takes 
Lj ~ 0.8 nH, Cj ~ 20 fF, one obtains u p ~ 2n x 40 GHz 
and a resonator frequency u>q ~ 2ir x 8.7 GHz. These values 
imply a phase drop \cf) r \ — 0.0835 such that the strength cou- 
pling of the QPC-resonator interaction at the resonance phase 
in Eq. (|6j reaches g/c^o ~ 0.21, well situated in the ultra- 
strong coupling regime of light-matter interaction. This value 
is not restrictive and one we may reach ratios g/ujQ > 0.2 for 
suitable parameters of the Josephson junction, and realize the 



regimes that have receive d a great theoretical and experimen- 
tal attention in last year d 23 l 25 tt 27 l 32 i Figures [4J; and |4}l show 
spectra for the USC coupling strength g/too ~ 0.21. The solid 
black lines in Fig. RH correspond to the spectrum calculated 
from Hamiltonian Eb, while the dashed green lines represent 
the Jaynes-Cummings model. The discrepancy between the 
two models is due to the Bloch-Siegert shift introduced by the 
counterrotating terms in the full Hamiltonian Q. This shift 
was experimentally observed in a circuit QED setup with a 
flux qurP2, where a ratio g/u>o = 0.1 has been reached. In 
this sense, we pave the way to study quantum optical proper- 
ties associated to the QPC-CWR coupling in circuit QED. 

Finally, it is noteworthy to mention that dissipative mecha- 
nisms affecting the coherence of the QPC have already been 
studied 3 5 . The theoretical estimated relaxation and dephasing 
rates are of the order of 10 8 Hz. However this value can be 
substantially decreased by proper control of the electromag- 
netic environment^. 

Conclusions. — We have presented a general frame for 
studying the quantum dynamics of a superconducting point 
contact galvanically coupled to a single-mode resonator. We 
have analyzed three possible configurations that could be ex- 
perimentally implemented and exhibit distinct cavity QED 
features as compared to semiclassical models. In particular, 
we have shown that the proposed QPC-CWR dynamics could 
reach the strong and ultrastrong coupling regimes of light- 
matter interaction, allowing the proposed setup to be consid- 
ered as a novel quantum device for circuit QED technology. 
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